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A non-band theory is developed to describe the superconductor-insulator (SI) transtition in 
resistively shunted, single Josephson junctions. The I — V characteristic is formulated by a Landauer- 
like formula and evaluated by the path-integral transfer-matrix method. The result is consistent 
with the recent experiments at around 80 mK. However, the insulator phase shrinks with decreasing 
temperature indicating that the single Josephson junction becomes all superconducting at absolute 
zero temperature, as long as dissipation is present. 
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The resistively shunted single Josephson junction(RS J) 
at T = is described classically by the second order dif- 
ferential equation for the phase difference between the 
two superconductors, 
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Here C is the capacitance, Rt is the normal state re- 
sistance of the junction, Iq is the critical current, and 
/ is the bias current. A mechanical analogy of a parti- 
cle moving in a tilted washboard potential with friction 
may be helpful. The I — V characteristic resulting from 
(|l|) was studied by Stewart 0. The capacitance brings 
about hysteresis. In the weak dissipation limit i?x — ► oo 
current is entirely carried by Cooper pairs at V = and 
up to the critical current. This is the dc Josephson effect. 
For I > Iq, the I — V curve switches to the Ohmic branch 
where current is entirely carried by electrons, V = RI. 
In the opposit limit i?x - * hysteresis disappears and 
supercurrent at V — is gradualy replaced by the nor- 
mal current with increasing voltage. The thermal and 
quantum fluctuation is expected to affect only the regime 
I/I < R q /R T = 9, where R q = h/Ae 2 = 6.45 KQ. is the 
quantum resistance. Now for a typical charging energy of 
IK, the critical current is Iq — 2eEj/% — A2Ej/Ec nA 
whereas the parameter ranges associated with the SI 
transition are g < 1, Ej/Ec 1. The range of current 
where fluctuation is important is thus up to a few nA. 
The focus of the present paper is precisely in this current 
regime. It should be noted that if g ^> 1, the bias current 
could be as large as the critical current while still staying 
in the regime where fluctuation is important. Indeed this 
case was studied theoretically and experimentally ||] 
focusing on the negative differential resistance as a signal 
of the Coulomb blockade of Cooper pair tunneling. 

A microscopic description is to take into account the 
quasi-particle tunneling but with zero superconducting 
gap (resistively shunted) as a source of dissipation. The 
path-integral expression of the partition function is 0] , 



Z = J V<t> exp(-S/h) 
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The single Josephson junction (JJ) has also been stud- 
ied by a phenomenological model || which is just to re- 
place the term cos ^( r )~^( T m (||) by 0(t)(/>(t') lead- 
ing, upon stationary phase approximation, to ([IJ) with 
retarded dissipation. 

The first question concerning the phase dynamics is 
whether it is coherent or incoherent. A band theory as- 
sumes a coherent motion of phase over many local min- 
ima of the washboard potential. In band picture, the 
voltage calculation is easy because it is proportional to 
the group velocity of the quasi-charge M . On the other 
hand, a non-band theory is appropriate for incoherent 
motion, and the problem is reduced to the calculation of 
the escape rate of particle out of a dissipative metastable 
well, the quantum Kramers rate. 

Watanabe and Haviland discussed their recent ex- 
periment on the Coulomb blockade at g 1 in light of 
the band theory. However, other experiments j?]]^] have 
not found the negative differential resistance in the SI 
transition regime g < 1. This indicates that the phase 
dynamics near the SI transition is dominated by inco- 
herent tunneling. Indeed a typical experimental setup 
is, E c = 1 K, Ej/E c = 0.1 - 20, T/E c = 0.1-1 
and g < 1. And most importantly, the minimum volt- 
age experimentally observable is of order 0.5fiV which 
translates to the temperature 10 mK. So when the en- 
ergy splitting due to tunneling exceeds this value, the 
corresponding state is experimentally recognized as insu- 
lating. However, 10 mK is much smaller than the poten- 
tial barrier Ej at the SI transition. Near the SI transition 
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at g = 0, therefore, the effect of many wells is expected 
unimportant. Switching on dissipation, the SI transition 
shifts to a smaller critical Ej/Ec, but then dissipation 
will suppress the effect of incresed tunneling strength. 
One thus expects that the phase dynamics essentially at 
a single well would determine the SI transition. 

In this paper we propose a non-band theory for the 
SI transition 0. Two methods have been employed to 
calculate the quantum Kramers rate along this line. One 
is the imaginary free energy method with an instanton 
technique [p|JlO|] . The other is the quntum transition state 
(QTS) theory long known in chemistry Here we 

develop a new approach combining a Landauer-like for- 
mula for the quantum Kramers rate |l3f| and the cluster 
transfer matrix (TM) method |0. The Landauer for- 
mula for electrical conductance has proven to be quite 
useful in mesoscopic systems. In the quantum Kramers 
rate problem, there is formally exact formula jy| anal- 
ogous to the Green-Kubo-Nakano formula for electrical 
conductivity. The corresponding Landauer-like formula 
is the quantum transient state (QTS) theory which tries 
to calculate the Kramers rate by summing up all the 
probabilities of finding particle at the top of the bar- 
rier which have positive, outgoing momentum. The QTS 
theory, however, has often been exercised with a semi- 
classical approximation |l^,[l6| and its drawbacks do not 
show up this way. As soon as one tries to exactly evaluate 
the path-integrals, however, one encounters difficulticss. 
First, the contribution to the escape rate is no longer 
local. The phase motion at different places than at the 
barrier top contributes. Moreover, the non-local contri- 
bution could become negative, as discussed in p3| . This 
undesirable nonlocality and negativity originates from 
the Weyle quantization procedure The original 

idea to collect all the outgoing, positive contribution at 
the barrier site, is thus spoiled. 

We have recently proposed a Landauer-like formula for 
the quantum Kramers rate. Applied to the voltage cal- 
culation here, it reads Q 

_Z_ = [i _ exp(-7r^/3/ e )] x ^{5(0-0o)P 2 /L) (5) 

where P = — 9^, </>o is the barrier top in the potential 
U(<f>), L is the system size, and < ... > is a thermal aver- 
age. The first factor in (5) takes into account the back- 
ward flux by detailed balance, whereas the second term 
is esentially square root of average of the kinetic energy, 
in mechanical analogy, at the barrier top with factor 1/2 
taking into account only the right-going contribution. In 
contrast, the QTS formula is obtained from (5) by replac- 
ing y/(6((/> - 4>o)P 2 /L) -> (Sty - <p )\P\)- The difference 
between the formula (5) and the QTS theory is thus a 
fluctuation which is precisely the undesirable non-local, 
possibly negative contributions in the QTS theory de- 
scribed in the above. After some manipulations, one can 
express the average in (||) as 



(6(0 - <f> )P 2 ) = Z/Z 



Z = J dq>W(cj), 0) 
in terms of the path integral 

W(<p, (j)') = I V4, exp(-S/h) 



(6) 
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where (j>(0) = 4> and 4>(J3ft,) = (/)'. 

The path-integral (g) with the action S given by (j^) 
can be evaluated precisely by the cluster TM method 
|jl4|| . In the present problem, however, the dimensionless 
conductance is at g < 1, and from the study of the sin- 
gle electron box which has a similar action as (||), this 
regime of g can be accurately handled by the 1-cluster 
TM method as follows. First divide the interval (0,/3K) 
into TV segments of size A = j3U/N . Denoting the dis- 
crete phase points as 4>ii i — 1, 2, iV + 1, we have 
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exp(-S/h) = J|x(0i,0 i+ i) x Rm 

8=1 

where the TM operator K is defined by 
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K = exp 
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with dimensionless parameters J = Ej/Ec, T 
k B T/E c , I = Ih/(2eJ) and 



Rm = exp 
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where the cumulant approximation was used which is pre- 
cise for g < 1. Second write 



r N+i r 



i S (4>n+i 0') (13) 



The first ^-function can be written in terms of a complete 
set {ip P } as 



(14) 



Now one can choose the complete set such that it satisfies 
the TM equation, 



d<f>K(<i>, <j>')ip p (<j)) = \ P ^ p (cj>) 



(15) 
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Repeatedly using the TM equation, we arrive at 

W(<j ) ^')=Rm-Y / ^r p WM^) (16) 



The correlation function (cos ^ J ~^' ) in ( |l2| ) can be calcu- 
lated likewise. 

There are three parameters T, g and Ej/Eq- In ad- 
dition the present method contains two more parame- 
ters A and Jr. The latter defines the phase space; 
4> a = — 7T — 4>o — tt I Jr> 4>b = T — 4>o + 7t/ Jr . For numeri- 
cal results obtained, we have checked that the results are 
affected a couple of % for the two different Jr = 2, 3 and 
for the two different A = 0.05,0.1. The results shown 
below are for Jr = 3 and A = 0.1. The "zero-bias" re- 
sistance reported in experiments may need to be treated 
carefully. In fact, the voltage measurment in a single J J is 
limited, unlike ID or 2D JJ arrays, to V m i n ~ O.bfiV M. 
So we have calculated the differential resistance R = 
with the constraint that the voltage is equal to O.BfiV 
(The calculation is repeated for V m i n = O.ObfiV to find 
no essential difference from the results below). At the 
same time, to be consistent with our theoretical consid- 
eration above, we have checked that the current is much 
smaller than the critical one, I -C Io- For results be- 
low, the current / is up to a few % of the critical one 
Io- Fig. 1 shows a temperature dependence of the resis- 
tance at V = 0.5^^ for g = 0.4. The small Ej/Ec case 
shows an insulator-like temperature dependence ^ < 0, 
while the large Ej / Ec case shows a superconductor- like 
behavior ^ > 0. Repeating the calculations for differ- 
ent g, one can draw a phase diagram. Fig. 2 is a phase 
diagram at T = 80 mK . The experimental results [0J|] 
are open circles (superconductor-like) and solid circles 
(insulator- like) at T ss 80 mK. Our result is denoted by 
open diamonds. The band theory result, thick solid line, 
is obtained by simply replacing the minimum measura- 
ment energy eV m i n by the thermal energy ksT |Q. The 
QTS theory result, triangles, is also plotted for compar- 
ison. A major disagreement between theory and experi- 
ment is for some data points near g = 2.8 ||. However a 
similar phase diagram experimentally found for the 2D JJ 
arrays with similar parameter ranges for <?, Ej/Ec an d T 
is bounded, E.,/E c < 0.5, and g < 0.5 (cf. Fig 3 in @). 
The above data near g = 2.8 is currently mysterious. 

There are two notable points in the obtained temper- 
ature dependence of resistance. First it becomes flat for 
large Ej/Ec irrespective of g. This is consistent with 
experimental findings 0(8). Theoretically, when the ra- 
tio Ej/Ec is sufficiently large, the quantum Kramers 
rate would become less dependent on T and g. Note a 
sharp contrast between the single JJ and the 2D JJ ar- 
rays |17|. The temperature dependence of resistance in 
the latter case showed a true phase transition behavior, 
namely resistance decreases or increases roughly expo- 
nentially with temperature. Interestingly, the tempera- 



ture dependence of resistance in ID JJ arrays also be- 
comes flat deep inside the superconducting phase p8| . 
The second notable point is that the SI phase boundary 
defined by vanishing dR/dT actually depends on temper- 
ature. In Fig. 3, we have plotted the SI phase boundary 
for T = 80, 60, 40, 20, 10 and 5 mK. Clearly, the insu- 
lator phase diminishes with decreasing temperature. We 
thus reach a conclusion that, as long as dissipation is 
present, the single JJ becomes all superconducting at ab- 
solute zero temperature. Note that, in contrast, the band 
theory predicts an essentially temperature independent 
phase boundary. It is interesting to see if our finding can 
be observed experimentally. 

Finally, related to the last point, we should argue that 
our finding is indeed consistent with the previous the- 
oretical works on a different RSJ model 10,01 and the 



subsequent experimental observations |19| of a universal 
T 2 -law for the escape rate k + of the zero- voltage state in 
the single JJ. That is, 



ln[k + (T) / k + (0)} cx T 2 for T -> 



(17) 



These experiments focused on the highly dissipative 
regime g> 1, and the current is very close to the criti- 
cal current, and therefore the escape rate k + is not quite 
proportional to voltage. In contrast, we focused on the 
regime g < 1 and calculated voltage since current is only 
up to a few % of the critical one. Nevertheless the present 
theory and the previous theories agree in that both k + 
and R(T) increase with temperature at very low T. In 
fact it was argued that the universal low temperature 
enhancement arizes from the temperature dependence of 
the heat-bath Jl^| . In the present case, it is due to the 
a(r) term in the action ([}]) and its temperature depen- 
dence originates from the Fermi distribution function in 
the electron Green's function |J. As T — > 0, this term 
dominates the temperature dependence and the voltage 
behaves as 



ln[V(T)/V(0)] cx gT 2 



(18) 



leading to ^Ji > for T — > 0. Note that, as is clear from 
( |l8| ) and seen in Fig. 3, such a temperature dependence 
disappears at g = 0. 
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FIG. 1. Temperature dependence of the " zero-bias" re- 
sistance at g = 0.4. From top to bottom, Ej/Ec=0.Q7, 
0.56, 0.84, 1.12, 1.4, 1.96, 2.52, 3.0 and 4.0. 
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FIG. 2. Phase diagram of shunted Josephson junc- 
tion at T — 80mK. The phase boundary lies between 
the insulator-like (solid circles) and superconductor-like 
(open circles) samples experimentally found in [7,8]. The 
thick line is the band theory (B-T), the triangle the QTS 
theory and the diamond is due to the present theory (L- 
like). 




FIG. 3. Temperature dependence of the 

superconductor-insulator phase boundary in the 
Ej/Ec — g plane. From right to left, T = 
80,60,40,20,10, and 5 mK. 
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